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摘 ”要 : SMAP(Soil Moisture Active Passive ,SMAP) 产 品 空间 分 辨 率 低 的 特征 限制 了 在 地 表 高 异 质 性 的 干旱 区 沙漠 


稀 玲 植被 区 的 适用 性 。 考 虑 到 干旱 区 沙漠 植被 区 特殊 的 环境 特征 ， 
一 化 植被 指数 (Normalized Difference Vegetation Index , NDVD)、 数 字 高 程 模型 (Digital Elevation Model, DEM) 等 传统 降 
尺度 因子 的 基础 上 ,增加 了 与 荒漠 地 表土 壤 水 分 关联 性 更 强 的 增强 型 修改 土壤 植被 指数 (Enhanced Modified Soil- 


在 地 表 温 度 (Land Surface Temperature ,LST)、 归 


Adjusted Vegetation Index,EMSAVD 与 比值 沙 地 亮度 指数 (Ratio Sand Brightness Index,RSBIT) 分 别 作为 反映 研究 区 植 


被 盖 度 和 裸 沙 分 布 状况 的 降 尺 度 因子 ,条 


=. 


用 随机 森林 算法 (Random Forest, RF), T ££ FARK HOOK AT BEN BEBS 


型 。 结 果 表 明 :(1) 由 相关 性 分 析 可 知 ,EMSAVI(rf= - 0.37, r w= — 0.34) RSBI (r 42 - 042, r2 — 0.25) 对 荒漠 土壤 
水 分 均 有 较 好 的 指示 作用 且 效 果 优 于 NDVI(r+=- 0.21,ra=0.08);(2) EMSAVI 和 NDVI 重 要 性 分 别 为 18.796 、 


13.2% ,EMSAVI 在 构建 降 尺度 模型 时 贡献 度 更 高 。(3) PA 


a 的 干 、 湿 季 干 旱 区 土壤 水 分 降 尺 度 模型 得 到 的 结果 与 


SMAP j^ f ÉY R4) 3A 8 T 0.916,0.910, RMSE 分 别 达到 了 0.0075 em? - em? 0.0063 em? cm”, 较 传统 模型 的 RMSE 
TEE IK T 0.0013 em? em, (4) 通过 计算 LBP(Local Binary Patterns) 的 差 值 (LBP.) 对 空间 一 致 性 评价 ,新 构建 降 尺 度 
模型 的 结果 (0.0585 ) 优 于 传统 降 尺 度 (0.0645)。 研 究 结果 将 短波 红外 波段 引入 到 植被 指数 建立 的 EMSAVI, 可 较 好 


地 应 用 于 干旱 区 沙漠 稀 琉 植被 区 土壤 水 分 降 太 度 研究 。 


关键 词 : 土壤 水 分 ;随机 森林 ; SMAP; MERE; 古 尔 班 通 古 特 沙漠 


土壤 水 分 (Soil Moisture, SM) 作为 全 球 水 循环 的 
重要 组 分 ,通过 蒸发 控制 着 大 气 圈 的 水 文学 循环 和 
气候 变化 ,是 控制 地 球 各 圈 层 内 部 水 能 量 交换 的 重 
要 参数 之 一 " ;同时 土壤 水 分 通过 生态 系统 水 胁迫 
来 降低 总 初级 生产 力 (Net Primary Productivity， 
NPP), 导致 植被 死亡 ,是 影响 陆地 矶 吸收 的 重要 
R 因此 ,准确 获取 土壤 水 分 信息 是 研究 陆地 生态 
系统 水 碳 循环 的 一 个 关键 问题 。 目 前 获取 土壤 水 
分 的 途径 有 直接 方法 和 间接 方法 。 直 接 方法 即 传 
统 土壤 水 分 观测 方法 ,主要 利用 探 针 或 重量 测量 法 
在 小 范围 内 测量 土壤 水 分 。 无 法 获取 大 范围 高 时 
相 土 壤 水 分 信息 。 相 较 于 计 重 采样 等 地 面 测量 的 
传统 方法 ,利用 卫星 遥感 探测 土壤 水 分 信息 的 手段 
具有 很 大 的 优势 。 其 中 微波 遥感 由 于 具有 穿 透 能 
力 强 与 土壤 水 分 具有 较 强 的 物理 联系 等 特点 ,而 
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成 为 监测 土壤 水 分 信息 的 主要 手段 之 一 *”*。 目 前 ， 
众多 卫星 平台 搭载 了 微波 传感器 且 基 于 微波 的 土 
壤 水 分 产品 也 达到 了 较 高 的 反 演 精度 ,如 主 被 动 观 
测 结合 的 SMAP 卫 星 .土壤 水 分 海洋 盐 度 SMOS E 
Æ „Sentinel 系列 主动 微波 卫星 以 及 AMSR2 .SMAP 
土壤 水 分 产品 等 。 基 于 微波 的 土壤 水 分 产品 往往 
时 间 分 辨 率 较 高 ,可 大 范围 .高 时 相 获 取 土 壤 水 分 
Fi AS, ,但 其 空间 分 辨 率 较 低 ( 数 千 米 到 数 十 千 米 )， 
不 能 很 好 的 表达 强 异 质 性 区 域 的 土壤 水 分 。 

被 动 微波 土壤 水 分 空间 降 尺 度 是 获取 高 时 空 
分 辩 率 土壤 水 分 产品 的 有 效 途 径 。 被 动 微波 降 尺 
度 方法 是 采用 模型 确定 低空 间 分 辨 率 微波 像 元 值 
与 高 空间 分 辨 率 地 表 信息 两 者 之 间 的 关系 ,并 借 此 
实现 空间 分 辨 率 的 提高 "。 有 研究 表明 , 降 尺 度 因 
子 往往 选择 NDVI, LST 以 及 数字 高 程 模型 DEM 作 
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为 降 尺 度 因 子 构建 经 验 降 尺度 模型 ** ,但 是 对 于 古 
尔 班 通 古 特 沙 漠 稀 玻 植 被 区 ,植被 指数 易 受 裸 沙 影 
响 ,严重 限制 了 对 土壤 水 分 关联 关系 的 表达 。 因 此 
对 于 干旱 区 沙漠 稀 足 植被 区 ,找到 一 些 与 土壤 水 分 
有 更 好 的 关联 关系 降 尺 度 因子 是 实现 降 尺 度 的 关 
键 环节 。 

本 文 以 古 尔 班 通 古 特 沙漠 为 研究 区 ,基于 随机 
森林 算法 ,在 优化 修改 型 土壤 调整 植被 指数 (Modi- 
fied Soil Moisture Vegetation Index , MSAVI) 基 础 上 提 
出 EMSAVI 指 数 ,构建 适用 于 沙漠 环境 特征 的 土壤 
水 分 降 尺 度 模 型 。 将 综合 考虑 稀疏 植被 和 裸 沙 的 
新 土壤 水 分 降 尺 度 模型 与 传统 降 尺 度 模型 进行 对 
比 ,探索 新 的 降 尺 度 模型 对 SMAP 土壤 水 分 产品 降 
尺度 的 合理 性 ,实现 高 时 空 分 辩 率 土壤 水 分 的 探 
测 。 以 期 为 干旱 区 土壤 水 分 降 尺 度 的 研究 提供 思 
路 ,为 荒漠 生态 系统 碳 汇 研 究 提供 重要 的 数据 
xd, 


1 材料 与 方法 


1.1 研究 区 概况 

古 尔 班 通 古 特 沙漠 (44”15'~46”5'N, 84°50" ~ 
91^20' E) f TF VERS ZR da Hb np p 33 28 39 du] Tj JS 4 7 
及 乌 伦 古河 以 南 地 区 (图 1)。 是 中 国 第 二 大 沙漠 ， 
面积 约 4.88x10' km? 7, Jt iln FERE ,是 新 
38 Jb op dtu EE AS ASE, EE ME et 70- 
150 mm, 受 冬季 融雪 水 人 渗 补 给 EFKA Zh 
给 和 蒸 散 发 的 共同 影响 ,春季 为 土壤 水 分 补给 期 
( 湿 季 ) ,夏季 为 土壤 水 分 耗损 期 ( 干 季 ) ,蒸发 量 约 
为 2000~2800 mm, 年 均 气 温 7.19 °C ,土壤 以 风沙 土 


46°N 


45°N 


44°N 


注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2020)4619 
号 ,对 底 图 边界 无 修改 。 下 同 。 
图 1 研究 区 概况 示意 图 
Fig.1 Schematic diagram of the study area 


为 主 ,容重 1.54~1.65 gem3035。 沙 漠 内 部 大 部 分 
为 固定 和 半 固 定 沙 丘 , 固 定 沙 丘 上 植被 覆盖 度 大 约 
3096-5096 , 半 固 定 沙 丘 大 约 达 15%~25%057。 融 雪 
水 和 春季 降水 作为 古 尔 班 通 古 特 沙漠 植被 生存 的 
主要 水 分 来 源 , 占 全 年 降水 的 65% 以 上 中 ;冬季 有 
稳定 积 雪 且 稳 定 积 雪 日 100~150 d, 积 雪 的 消融 下 渗 
为 春季 短命 植物 的 生长 提供 了 有 利 的 条 件 , 受 此 影 
响 春 夏季 短命 植物 生长 旺 威 **。 

1.2 数据 来 源 

1.2.1 土壤 水 分 数据 SMAP 土壤 水 分 数据 采用 
SMAP 增强 型 全 球 每 日 9 km 土壤 水 分 产品 (SMAP 
Enhanced L3 Radiometer Global Daily 9 km EASE- 
Grid Soil Moisture, V004, F Xf ff Jj SMAP L3) , 反 
Ilt 0-5 cm 土壤 深度 的 平均 含水 量 ,单位 为 cm em ^, 
上 午 ( 降 轨 ) 和 下 午 ( 升 轨 ) 的 星 下 点 过 境 时间 分 别 
为 当地 时 间 早 上 06:00 和 下 午 18:00。 由 于 近 地 表 
土 层 和 植被 覆盖 的 热平衡 和 均 质 大 气 环境 在 上 午 
的 时 候 更 接近 地 表 真 实情 况 ” ,选择 该 产品 的 降 轨 
数据 作为 粗 分 辨 土壤 水 分 。 通 过 降水 量变 化 和 研 
究 区 已 有 的 研究 结论 ”, 并 考虑 SMAP L3 数据 覆盖 
完整 情况 ,分 别 选取 2018 年 湿 季 . 干 季 中 4 月 15 日 
和 8 月 24 日 的 土壤 水 分 作为 检验 新 旧 降 尺度 模型 
估算 土壤 水 分 的 验证 数据 ,探索 方法 的 可 行 性 。 
1.2.2 逐 感 数据 MODIS 数 据 以 MODO9A1 的 8d 
合成 反射 率 产 品 为 基础 ,计算 NDVI、RSBI、EMSAVI 
数据 ;地 表 温 度数 据 选用 MOD11A1 日 地 表 温 度 产 
im (https://ladsweb.modaps.eosdis.nasa.gov ) o 

1.2.3 x75 44% SRTM(Shuttle Radar Topography 
Mission), 由 美国 太空 总 署 (NASA) 和 国防 部 国家 测绘 
局 (NIMA) 联 合 测量 , HES HE HE 60° N~56°S 的 高 分 辩 
率 地 表 高 程 数 据 。 研 究 使 用 的 是 SRTM DEM 30 m 
分 状 率 高 程 数据 (https://earthexplorer.usgs.gov/)。 

1.3 研究 方法 

1.3.1 随机 森林 算法 众多 人 研究 已 经 讨论 了 集成 学 
习 方 法 (如 随机 和 森林) 来 解决 遥感 中 的 分 类 和 非 线性 
[n] I [a] 7 | BAL AR AK CRE) H Breiman ”提出 ， 
是 一 种 用 途 广 泛 ,高度 灵 活 的 机 器 学 习 算 法 ,通过 
bootstrap 重 采 样 技术 从 原始 样本 中 抽取 并 生成 训练 
PEASE AB ,然后 根据 训练 样本 集 生 成 多 个 决策 树 ， 
最 终 通过 决策 树 的 投票 分 数 确 定 结果 站 ” 。 根 据 算 
法 的 自 适 应 性 质 ,RF 在 确定 自 变量 和 因 变 量 之 间 复 
杂 的 非 线性 关系 方面 具有 优势 ” ,与 其 他 机 器 学 
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习 方法 相 比 ,随机 和 森林 抗 过 拟 合 能 力 强 ,异常 值 和 
噪声 的 容忍 度 高 ”。 此 外 ,模型 内 置 了 特征 重要 性 
评价 机 制 ,拥有 较 强 的 模型 可 解释 性 ,变量 重要 性 
评分 (Variable Importance Measures, VIM) 用 以 计算 
每 个 特征 在 随机 森林 中 的 每 棵 树 中 所 做 出 的 贡献 ， 
通常 选择 基尼 系数 (Gini index) 或 者 袋 外 数据 误差 
(Error of out-of-bag, Eoo) 进行 衡量 ,评价 公式 如 下 : 
上 oop = wo’) = Y] (1) 
式 中 :N 为 00B 样 本 个 数 ; Y QC) 为 给 定 样 本 X, 时 ， 
RF 模型 的 预测 数据 ;7; 为 实际 数据 。 
首先 计算 每 棵 决策 树 Eoo ,然后 对 特征 M 随 
机 加 入 噪声 ,并 再 次 计算 袋 外 数据 误差 Eoos, ,由 此 
得 到 特征 M 的 重要 性 ( VIM, ) 计 算式 : 
VIM, = HS Eson = Eige) (2) 
式 中 :为 随机 森林 中 的 树 的 数量 。 
若 给 某 个 特征 随机 加 入 噪声 之 后 , 袋 外 数据 误 
差 大 幅度 提高 , 则 说 明 这 个 特征 对 于 数据 的 预测 结 
影响 很 大 , 即 此 特征 重要 性 较 高 。 因 此 ,本文 基 
于 随机 和 森林 算法 整合 多 元 参数 来 反 演 土壤 水 分 。 
1.3.2 土壤 水 分 降 尺 度 模 型 构建 ”目前 ,随机 森林 
算法 的 土壤 水 分 降 尺 度 模 型 ,往往 是 基于 土壤 水 分 
与 植被 盖 度 、 地 表 温 度 三 者 之 间 关 系 的 理论 基础 ”。 
与 土壤 水 分 有 密切 联系 的 辅助 数据 作为 降 尺 度 因 
子 , 如 地 表 温 度 .植被 指数 .地 形 因子 等 ,其 中 NDVI 
和 地 表 温 度 被 广泛 使 用 。 然 而 ,相关 研究 表明 , 土 
壤 水 分 不 仅 受到 植被 盖 度 影响 ,而 且 受到 地 表 善 度 
类 型 和 土壤 属性 的 影响 ”。 故 在 传统 降 尺 度 模型 
中 降 尺 度 因子 的 基础 上 ,考虑 稀 叭 植被 和 裸 沙 的 环 
境 特征 ,增加 了 与 荒漠 土壤 水 分 关联 性 更 高 的 EM- 
SAVI 与 RSBI 分 别 作为 反映 研究 区 植被 盖 度 和 裸 沙 
分 布 状况 的 降 尺 度 因 子 ,构建 干旱 区 土壤 水 分 降 尺 
度 模 型 。 
(1) 增强 型 修改 土壤 植被 指数 (EMSAVD) 
植被 指数 利用 光谱 信息 实现 对 植被 状态 信息 
的 表达 ,如 NDVI 可 定量 地 评价 植被 盖 度 .生长 活力 
及 生物 量 等 。 但 NDVI 对 土壤 和 植被 的 敏感 度 不 
同 ,在 植被 密度 低 的 地 区 ,NDVI 受 土壤 背景 影响 较 
大 ,会 忽略 部 分 植被 信息 ;在 植被 密度 高 的 地 区 , 
NDVI 则 会 出 现 光 谱 饱 和 现象 。 而 改进 的 植被 指数 
(MSAVD 整 合 了 土壤 背景 状况 和 大 气 调整 因子 ,使 
用 MSAVI 能 使 植被 指数 固有 的 背景 “噪声 " 降 到 最 
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小 ,便于 消除 植被 冠 层 信息 及 土壤 背景 的 影响 ” 。 
2BNm Tl- JOB vis + 17 "E SB -B,) 
2 


MSAVI = 


(3) 
SUF Bw ZEE BEN c8] 2 ; Br AJ £L EE BS c 
射 率 。 有 研究 表明 ,只 有 2 个 波段 组 成 的 植被 指数 
因 有 显著 的 相关 性 ,而 导致 构建 的 模型 不 稳定 2 ， 
故 提 出 EMSAVI。 
2B. + 1- [(2Bym + 1 - 8B — Be) + Bas 
2+ Bean 


EMSAVI= 


(4) 
式 中 :Brn 为 近 红 外 波段 的 反射 率 ;Br 为 红 波 段 的 反 
射 率 ;Bswir 为 短波 红外 波段 的 反射 率 。 有 研究 发 现 ， 
短波 红外 波段 对 土壤 水 分 变化 具有 重要 指示 作用 
且 对 土壤 水 分 反 演 精度 提升 有 一 定价 值 ” 且 与 可 
见 光 和 近 红 外 相 比 ,Bswm 对 背景 的 敏感 度 较 小 二 , 故 
在 MSAVI 基 础 上 引入 对 土壤 水 分 敏感 度 更 高 的 短 
波 红外 波段 ,得 到 EMSAVI, 减 少 植被 冠 层 及 沙漠 对 
土壤 水 分 探测 的 影响 ,最 终 提高 植被 指数 与 土壤 水 
分 的 相关 性 。 

(2) 比值 沙 地 亮度 指数 (RSBI) 

RSBI 指 数 综合 红 波段 .短波 红外 与 绿 波段 3 个 
波段 的 光谱 信号 ,增强 了 光谱 信息 ,该 指数 通过 计 
算 沙 地 反射 率 最 高 与 最 低 波 段 间 的 比值 ,进一步 扩 
大 波段 间 差异 范围 ,从 而 在 影像 上 提高 了 沙 地 与 非 
沙 地 识别 的 灵敏 度 ”。 

RSBI = Bswm +B, 
B 


式 中 :Bswr 表 示 短 波 红外 波段 的 反射 率 ;B 表 示 红 波 
段 的 反射 率 ;Be 表示 绿 波段 的 反射 率 。 

利用 随机 森林 采用 非 线性 函数 本 算法 建立 EM- 
SAVI, NDVI, RSBI, DEM 及 LST 与 土壤 表层 水 分 
(SW) 之 间 的 相互 关系 : 

SM=F,(EMSAVI, NDVI, RSBI, DEM, LST) (6) 

基于 SMAP 被 动 微 波 土壤 水 分 的 干旱 区 表层 土 
壤 水 分 降 尺 度 技 术 路 线 如 图 2 所 示 。 

降 尺 度 模 型 构建 步 又 如 下 : 

CL) 1 km4% A EMSAVI, NDVI, RSBI, LST 
FI DEM 数据 重 采样 至 9 km 分 辨 率 , 与 9 km 空间 分 
辩 率 的 SMAP 土壤 水 分 数据 的 像 元 进行 匹配 。 

(2) 将 9 km 空间 分 辩 率 的 SMAP 土壤 水 分 数据 
作为 目标 , 升 尺度 后 得 到 的 EMASVI , NDVI, RSBI, 
LST 5j DEM 数据 作为 降 尺度 因子 ,组 成 2 种 组 合 使 


(5) 
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训练 模型 


图 2 技术 路 线 图 


Fig. 2 Flowchart of the downscaling procedure in this study 


用 随机 森林 进行 训练 ,在 9 km 空间 尺度 上 构建 降 尺 
度 因子 与 土壤 水 分 之 间 的 映射 关系 ,利用 土壤 水 分 
数据 对 新 旧 降 尺度 模型 进行 精度 验证 并 对 比 , 检 验 
新 构建 模型 在 干旱 区 荒漠 的 适用 性 并 通过 一 些 评 
价 指标 对 精度 验证 结果 进行 评价 。 

(3) 将 高 分 辨 率 降 扩 度 因子 重 采 样 至 1 km 空 
间 分 辨 率 ,并 应 用 于 降 尺 度 模 型 中 ,得 到 1 km 空间 
分 辨 率 土壤 水 分 数据 。 
1.3.3 评价 指标 ”本 人 研究 选取 一 些 典 型 的 评价 指标 
对 基于 随机 和 森林 的 新 旧 土 壤 水 分 降 尺 度 模型 进行 
精度 评价 ,指标 包括 :决定 系数 R(Coefficient of De- 
termination ) 、 均 方 根 误差 (Root Mean Square Error, 
RMSE) 、 局 部 二 值 模式 (Local Binary Pattern ,LBP) 和 和 
BEAK AR 4H 2& AK BL (Pearson Correlation Coefficient, 
r)。 计 算 公 式 如 下 : 


Sm, = SM, 
i=l 


gn (7) 
S (SM, - SM) 
=1 
N 
S\(SM,- SM, 

= i=1 
RMSE = |=", (8) 
LBP = decimakX, X,, * * *, X) (9) 

SM, - SM,XSM, - SM, 

js > R aX L p) (10) 


JY (SM, - SM,’ (SM, - SM,y 
式 中 : SM, 为 土壤 水 分 预测 值 ; SM, 为 土壤 水 分 预 


测 平均 值 ;SM 为 原始 的 SMAP L3 土壤 水 分 数据 ; 
SM, 为 SMAP L3 土壤 水 分 平均 值 ;N 为 样本 数量 ;XX 
为 3x3 移动 窗口 中 中 心 像 元 周围 像 元 的 值 。 公 式 
如 下 : 


LBP, = S|SH, -SY (11) 
"Ex 


StH SHAN SY 4) 3g 9 km 降 尺 度数 据 和 9 km 土壤 
水 分 产品 的 LBP 值 ; LBP. 为 降 尺 度数 据 与 原始 产 
品 像 元 差 值 的 绝对 值 之 和 。 


2 结果 与 分 析 


24 降 尺 度 因子 与 土壤 含水 量 的 相关 性 分 析 

对 5 个 降 尺 度 因子 与 表层 (0~5 cm) 土 壤 水 分 计 
算 斯 皮尔 曼 相关 系数 及 显著 性 检验 ,进而 分 析 干 、 
湿 季 土壤 含水 量 与 不 同 降 尺度 因子 之 间 的 相关 性 。 
如 表 1 所 示 , 湿 季 士 壤 含水 量 与 各 降 尺 度 因子 相关 性 
为 :EMSAVI(-0.34) >LST(-0.31) >RSBI ( - 0.25) > 
NDVI(0.08)>DEM (0.07) ,而 传统 模型 中 降 尺 度 因 子 
之 一 的 NDVI 相 比 增加 的 光谱 指数 表现 出 较 弱 相关 
性 (r=0.08) , 低 于 EMSAVI(-0.34) 与 RSBI(-0.25 ) , 
日 与 土壤 水 分 间 无 显著 相关 ; 干 季 土壤 含水 量 与 各 
降 尺 度 因子 相关 性 为 :DEM(-0.59)>LST(-0.47)> 
RSBI(-0.42)>EMSAVI(-0.37)>NDVI(-0.21)。 F., 
湿 季 的 降 尺 度 因子 相关 系数 与 改进 后 降 尺 度 因 子 
相 比 ,LST 具 有 较 强 的 相关 性 ,产生 该 现象 的 主要 原 
因 是 地 表 温 度 对 水 分 蒸发 起 决定 性 作用 ,通常 情况 
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表 1 降 尺 度 因子 与 土壤 含水 量 相关 系数 
Tab. 1 Correlation coefficient between environmental 


variables and soil water content 


降 尺 度 因子 湿 季 干 季 
传统 降 尺 度 因子 NDVI 0.08 -0.21" 
DEM 0.07 -0.59" 
新 增 降 尺 度 因子 LST -0.31" -0.47" 
RSBI -025" -042" 
EMSAVI -0.347 -0.37" 


注 :** 表 示 在 P<0.01 水 平 显 著 。 


下 地 表 高 温 区 域 (或 神 沙 区 域 ) 土 壤 水 分 含量 较 低 ; 
DEM 在 干 季 时 表现 较 强 的 相关 性 ,原因 可 能 是 该 降 
尺度 因子 是 通过 对 气温 和 植被 的 影响 来 间接 影响 
土壤 水 分 含量 , 干 季 时 降水 较 少 继而 导致 DEM 对 表 
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干 、 湿 季 降 尺度 因子 重要 性 平均 值 /% 


LST | DEM EMSAVI NDVI RSBI 
降 尺 度 因 子 


图 3 和 干 、 湿 季 降 尺度 因子 重要 性 得 分 平均 值 


Fig. 3 Average of downscaling factor importance scores for 


two periods 


层 土壤 水 分 的 影响 大 幅 提 高 ;NDVI 在 干旱 区 沙 演 
稀 玻 植被 区 对 土壤 水 分 的 指示 作用 较 小 ,表明 传统 
降 尺 度 模 型 在 荒漠 地 区 适用 性 较 弱 。 

综 上 所 述 ,本 模型 新 增 的 EMSAVI(Cria=-0.34、 
r 42-0.37) . RSBI(ra=-0.25 +=-0.42) 均 高 于 传统 
模型 中 NDVI 与 土壤 水 分 的 相关 系数 (ria=-0.08、 
r= -0.21) ,表明 添加 的 EMSAVI 和 RSBI 对 荒漠 土 
壤 水 分 有 较 好 指示 作用 且 效 果 优 于 NDVI, 选 择 的 
降 尺 度 因 子 具 有 合理 性 。 不 同 降 尺 度 因子 与 表层 
土壤 水 分 的 相关 性 表现 出 不 同 的 稳定 性 ,EMSAVI、 
RSBI 和 LST 在 干 湿 季 中 与 表层 土壤 水 分 的 相关 系 
数 结果 中 具有 较 强 的 稳定 性 ,DEM 和 NDVI 则 受 干 
湿 季 影响 较 大 ,稳定 性 较 差 。 
22 降 尺 度 因子 重要 性 分 析 

在 非 线性 关系 构建 过 程 中 ,随机 森林 算法 能 
根据 输入 因子 即 降 尺度 因子 建立 模型 ,并 提供 每 个 
降 尺 度 因子 的 重要 性 。 重 要 性 武大 ,说 明 该 变量 对 
模型 的 影响 越 大 。 由 图 3 可 知 , 干 、 湿 季 降 尺度 因子 
重要 性 总 体 均 值 表现 为 :DEM (39.3% ) >LST 
(18.9% ) >EMSAVI (18.7% ) >NDVI (13.2% ) >RSBI 
(10.096) 

降 尺 度 因 子 与 土壤 水 分 相关 性 统计 结果 与 随 
机 森林 中 降 尺 度 因子 的 重要 性 结果 有 所 差异 。 干 
季 的 DEM(r=-0.59) 略 高 于 LST(r=-0.47) ,但 在 相同 
时 期 DEM 重要 性 相 比 LST 增 加 了 20.4% ;相关 性 与 
重要 性 分 析 结 果 均 表明 新 增 的 EMSAVI(r= -0.35, 
重要 性 =18.7% ) 对 新 构建 干旱 区 土壤 水 分 降 尺 度 模 
型 的 贡献 高 于 NDVI(r= -0.15 ,重要 性 =13.2% ) ;RS- 
BI 的 贡献 度 略 低 于 NDVI, 但 对 模型 的 贡献 度 依然 


较 高 (重要 性 =10.0% ) 。 结 果 说 明 ,相关 性 较 低 的 
NDVI 也 能 为 干旱 区 土壤 水 分 降 尺 度 模 型 的 构建 提 
供 一 些 信息 ,但 是 改进 后 得 到 的 EMSAVI 更 适合 作 
为 干旱 区 土壤 水 分 降 尺 度 因 子 构建 模型 。 
2.3 SMAP 土 壤 水 分 降 尺 度 新 旧 模 型 对 比 

由 于 人 研究 区 空间 异 质 性 明显 ,土壤 水 分 在 不 同 
季节 相差 显著 ,因此 分 别 选取 土壤 水 分 相对 差异 较 
大 时 间 节 点 中 云 量 较 少 的 降 尺 度数 据 与 SMAP L3 
精度 验证 后 对 比 ,以 此 检验 新 旧 模型 在 研究 区 的 适 
用 性 。 

图 4 为 2018 年 湿 季 和 2018 年 干 季 分 别 利用 传 
统 与 新 构建 的 降 尺 度 模 型 得 到 土壤 水 分 预测 值 与 
SMAP L3 的 散 点 图 。 由 图 4 可 知 , 干 、 湿 季 预 测 值 与 
产品 值 均 具有 和 较 高 的 决定 系数 (R>0.88) JH OR HE 
(r>0.93) 以 及 较 低 的 误差 (RMSE<0.009 cm em?) o 
在 湿 季 相 比 于 传统 降 尺 度 模型 ,本 文 构建 的 改进 降 
尺度 模型 精度 有 所 提升 ,提高 了 17.10% ;在 干 季 相 
比 传 统 降 尺度 模型 ,本文 构 建 的 改进 降 尺 度 模型 精 
EEFE Y 14.77%. RADHA , SIRF RE Bt 
植被 和 裸 沙 的 环境 特征 所 构建 的 改进 模型 预测 效 
果 , 整 体 相 较 于 传统 模型 略微 提高 了 精度 ,表现 出 
更 高 的 决定 系数 、 相 关 性 以 及 更 低 的 误差 ,更 适用 
于 干旱 区 的 土壤 水 分 降 尺 度 研 究 。 
2.4 降 尺度 结果 空间 一 致 性 评价 

SMAP 土壤 水 分 产品 与 基于 RF 的 回归 方法 建 
立 的 干旱 区 环境 因子 之 间 的 定量 关系 模型 性 能 
好 ,可 将 高 分 辩 率 变量 应 用 于 该 模型 ,从 而 得 到 高 
分 辩 率 的 表层 土壤 水 分 数据 。 本 文 对 1 km 空间 分 
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图 4 不 同时 


0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 
SMAP L3/(cm^cm?) 


期 土壤 水 分 降 太 度数 据 与 SMAP 土壤 水 分 散 点 图 


Fig.4 Soil water downscaling data in different periods and SMAP soil water scatter plot 


辩 率 表层 土壤 水 分 降 尺 度数 据 与 9 km 2 [R] 4) ES 
原始 SMAP 土 壤 水 分 数据 进行 了 空间 一 致 性 评价 及 
可 视 化 比较 (图 5)。 

评估 降 尺 度 结果 表征 SMAP 单 个 像 元 土壤 水 分 
空间 变化 的 能 力 是 十 分 重要 的 ,空间 精度 可 以 通过 
计算 9 km 降 尺 度数 据 和 土壤 水 分 产品 之 间 的 空间 

寺 征 (如 纹理 ) 的 差异 来 量化 。 因 此 ,本 文通 过 计算 
新 旧 模 型 9 km 降 尺 度数 据 与 土壤 水 分 产品 之 间 局 
部 二 值 模式 (Local Binary Patterns, LBP) 的 差 值 (LB- 
Pe) 评 价 新 旧 模 型 数据 与 SMAP L3 的 空间 一 致 性 ( 表 
2)。 总 体 结果 表明 ,干旱 区 土壤 水 分 降 尺 度 模 型 
(LBPc,=0.0585 ) 性 能 优 于 传统 降 尺 度 模 型 (LBP.。= 
0.0645) 。 

人 研究 结果 表明 ,改进 后 的 降 尺 度 模型 能 够 将 
SMAP 被 动 微波 土壤 水 分 产品 的 空间 分 辨 紊 提高 (9 
km 到 1 km), 改 善 了 原始 产品 的 空间 细节 ,同时 保证 
降 尺 度 前 后 数据 在 大 多 数 区 域内 土壤 水 分 的 空间 
格局 一 致 。 湿 季 时 ,高 含量 土壤 水 分 主要 分 布 于 东 


部 区 域 ; 干 季 时 ,高 含量 土壤 水 分 主要 分 布 于 西部 
和 西南 部 ,而 低 含量 土壤 水 分 出 现在 中 部 和 东南 部 
Kiko SMAP 土壤 水 分 产品 的 数值 范围 与 得 到 的 土 
壤 水 分 数据 基本 保持 一 致 (0.02~0.13 cm * em?) , 表 
明 与 其 他 基于 机 融 学 习 的 降 尺 度 算法 相 比 ,RF 能 够 
保持 粗 分 状 率 土壤 水 分 的 动态 范围 。 


3 讨论 


前 人 研究 主要 通过 NDVI 表 达 植 被 状态 信息 ， 
将 其 作为 一 个 重要 的 降 尺 度 因 子 构建 土壤 水 分 降 
尺度 模型 “”。 本 人 研究 区 为 古 尔 班 通 古 特 沙漠 ,植被 
稀 玖 ,传统 土壤 水 分 降 尺 度 模型 中 NDVI 因 子 易 受 
植被 密度 、 裸 沙 的 影响 ,在 干旱 区 土壤 水 分 降 尺 度 
研究 表现 出 极 大 的 局 限 性 ,因此 本 研究 基于 随机 森 
林 算 法 ,在 以 往 研 究 中 广泛 选择 LSTNDVI.DEM 作 
为 降 尺 度 因子 的 基础 上 ,同时 考虑 本 研究 区 的 土 
地 覆 被 类 型 的 差异 ,新 增 EMASAVI、RSBI 以 提高 干 


部 和 东南 部 ,而 低 含量 土壤 水 分 出 现在 中 部 和 西南 


早 区 中 植被 与 裸 沙 的 识别 精度 ,并 反映 裸 沙 对 表面 
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(a) 2018-04-15 SMAP L3 


(c) 2018-04-15 改进 模型 降 尺度 结果 


ChinaXiv 合 作 期 刊 


43 PEE ESE: TH AS BEES TRAE SMAP 土壤 水 分 产品 降 尺 度 分 析 589 


(b) 2018-08-24 SMAP L3 


(d) 2018-08-24 改进 模型 降 尺 度 结果 


(e) 2018-04-15 传统 模型 降 尺度 结果 


SM/(cm?-cm?) 
E Ee 


0.02 


(f) 2018-08-24 传统 模型 降 尺 度 结果 


0 62.5 km 
[m 


0.13 


图 5 古 尔 班 通 古 特 沙漠 干 di SMAP L3 和 相应 的 土壤 表层 水 分 降 尺 度数 据 


Fig. 5 SMAP L3 and corresponding soil surface water downscaling data for two periods in the Gurbantunggut Desert 


表 2 空间 一 致 性 评价 结果 


Tab. 2 Spatial consistency evaluation results 


湿 季 TÆ 
LBP.., 0.065 0.067 
LBP... 0.057 0.060 


IE: LBP..,LBP- 分 别 为 传统 模型 和 新 构建 模型 空间 一 致 性 评价 结 
果 。LBP. 值 为 0 时 ,表示 完美 的 融合 图 像 ; 越 接近 于 0 的 值 表示 空间 
一 致 性 越 好 。 


土壤 水 分 的 影响 ,建立 了 适用 于 干旱 区 的 表层 土壤 
水 分 的 降 尺 度 模 型 。 本 文 构建 的 干旱 区 地 表土 壤 
水 分 降 尺 度 相 较 以 往 传统 模型 的 误差 有 上 略微 的 降 
低 , 但 是 干 湿 季 精 度 均 提升 了 10% 以 上 。 本 文 在 干 
旱 区 地 表土 壤 水 分 降 尺 度 研究 中 的 降 尺 度 因 子 选 
择 及 地 表土 壤 水 分 降 尺 度 模型 的 构建 ,可 以 作为 后 
续 对 干旱 区 的 土壤 水 分 降 尺 度 人 研究 的 参考 ,如 : (1) 
本 文 在 MSAVI 中 引入 波长 范围 为 1.1~2.5 km, 能 有 
效 探测 植被 和 土壤 水 分 变化 ,对 背景 的 敏感 度 较 小 
且 包 含 更 多 信息 的 SWIR 波段 中 ,因此 引入 SWIR 
波段 可 以 大 幅 提高 植被 指数 与 土壤 水 分 之 间 的 关 
联 度 ,这 与 姚 云 军 等 2 的 研究 一 致 ;2) 因为 不 同 波 


段 或 光谱 指数 有 较 强 相关 性 汪 ,因此 本 文 在 MSAVI 
的 基础 上 加 入 SWIR 波段 ,可 以 最 大 程度 的 降低 波 
段 之 间 的 匈 余 信息 进而 提高 光谱 指数 的 指示 效果 ， 
这 与 蔡 亮 红 等 的 研究 一 致 ,日 在 不 同时 段 也 具有 
更 好 的 稳定 性 。(3) 加 入 的 RSBI 结 合 了 红外 、 近 红 
外 及 绿 波段 的 光谱 信息 ,能够 较 好 表示 沙 地 与 非 沙 
地 的 分 布 特征 3”, 故 构建 的 干旱 区 荒漠 土壤 水 分 降 
尺度 模型 能 够 将 SMAP 土壤 水 分 数据 由 9 km 空间 
分 辨 京 提升 至 1 km, 一定 程度 上 改善 了 SMAP 土壤 
水 分 数据 缺失 的 状况 ,最终 得 到 全 天 候 .高 时 空 分 
辩 率 的 土壤 水 分 数据 ,拓宽 随机 森林 算法 应 用 领 
域 ,为 数据 稀缺 区 表层 土壤 水 分 获取 提供 有 效 方 
法 ,说 明 本 文 构建 的 降 尺 度 模型 具有 较 大 潜力 ,可 
为 后 续 土 壤 水 分 含量 预测 提供 参考 。 

本 文 表明 基于 随机 森林 建立 的 土壤 水 分 降 尺 
度 模 型 能 够 建立 与 5 个 降 尺 度 因子 的 非 线性 关系 ， 
并 且 构 建 模型 具有 较 强 的 鲁 棒 性 ,这 与 im 等 % 的 研 
究 结 果 相 一 致 ;研究 发 现 地 表 温 度 对 土壤 水 分 相关 
性 及 对 降 尺 度 模型 的 贡献 率 均 较 高 ,是 影响 地 表土 
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培 水 分 的 关键 降 太 度 因 子 , 这 与 Zhao 等 ”的 研究 结  (〈-0.21)。 


果 一 致 。 在 以 往 研 究 中 ,DEM 与 地 表 水 分 具有 较 低 
的 相关 性 及 重要 性 ,但 在 本 研究 结果 中 海拔 表现 出 
较 高 的 重要 性 ,分 析 其 原因 可 能 是 沙漠 边缘 海拔 相 
对 较 高 的 绿洲 表层 水 分 数值 较 高 ,故而 提高 了 海 拨 
的 重要 性 。 在 本 文中 NDVI 虽 然 与 地 表土 壤 水 分 相 
关 性 不 高 ,但 是 在 构建 的 地 表土 壤 水 分 降 尺 度 模 型 
中 表现 出 一 定 的 重要 性 ,原因 可 能 是 由 于 研究 区 面 
IL Rig) ,地表 覆 被 差异 明显 ,NDVI 能 够 弥补 EMSA- 
VI 在 描述 特定 区 域 植物 生长 状况 方面 的 不 足 。 基 
于 随机 森林 的 降 尺 度 模型 相 较 于 传统 的 经 验 模型 
构建 更 加 便捷 ,并 且 能 够 提供 构建 模型 时 各 降 尺 度 
因子 的 重要 性 ,便于 对 环境 因子 与 土壤 水 分 关系 的 
逆向 推理 ,大 大 提高 了 获取 大 范围 .高 时 相 地 表土 
壤 水 分 的 可 行 性 ,可 为 干旱 区 表层 土壤 水 分 降 尺 度 
模型 研究 提供 参考 。 借 鉴 了 Zhu 等 “用 时 空 融合 方 
法 的 评价 指标 来 衡量 得 到 的 高 分 辩 率 水 分 的 精度 ， 
结果 表明 改进 的 降 尺 度 模 型 较 传统 模型 具有 更 高 
的 空间 一 致 性 ,更 适用 于 干旱 区 土壤 水 分 降 尺 度 的 
研究 ,后 续 可 以 增加 更 多 与 地 表土 壤 水 分 相关 的 数 
据 作为 模型 的 降 尺 度 因子 ,如 后 向 散射 系数 和 地 表 
反照 率 , 同 时 可 以 继续 提高 降 尺 度 结果 的 空间 分 辨 
率 , 以 期 更 好 的 减少 误差 和 克服 空间 异 质 性 。 


4 结论 


=A 


SMAP 是 可 以 表征 土壤 水 分 的 高 时 间 分 辩 率 的 
栅 格 数据 产品 ,但 较 低 的 空间 分 辩 率 限制 了 在 地 表 
高 异 质 性 的 干旱 区 沙漠 稀 玻 植被 区 的 适用 性 。 本 
文 利 用 随机 森林 算法 ,以 LST, NDVI, RSBI, DEM, 
EMSAVI 为 辅助 数据 ,探索 新 的 SMAP 降 尺度 模型 在 
干旱 区 沙漠 稀 玲 植被 区 的 合理 性 ,得 到 以 下 结论 : 

(1) 通过 相关 性 分 析 可 知 ,改进 的 EMSAVI 
( -0.35 ) 和 新 增 的 RSBI(-0.34) 均 高 于 NDVI 与 土壤 
水 分 的 相关 性 (-0.15) , 且 提 升 明 显 。 新 增 的 植被 指 
数 和 土壤 指数 对 匠 漠 土壤 水 分 有 更 好 的 指示 作用 ， 
效果 优 于 NDVI。 不 同 降 尺度 因子 与 地 表 水 分 表现 
出 不 同 的 相关 程度 ,在 不 同时 期 相关 程度 也 存在 差 
异 。 各 环境 因子 与 土壤 水 分 的 相关 性 表现 为 湿 季 : 
EMSAVI( -0.34)>LST(-0.31)>RSBI(-0.25 )>NDVI 
(0.08) >DEM (0.07) ; 干 季 表 现 为 :DEM(-0.59 )> 
LST(-0.47)>RSBI(-0.42 )>EMSAVI(-0.37)>NDVI 


(2) 降 尺 度 因子 相关 性 统计 结果 与 重要 性 结 
有 所 差异 。 干 季 DEM 相关 性 (r=-0.59) 略 高 于 LST 
(r=-0.47) ,但 同时 期 DEM 重要 性 相 比 LST 增 加 了 
20.4% ;在 环境 因子 相关 性 和 重要 性 结果 中 ,EMSA- 
VI 的 统计 结果 均 要 高 于 NDVI; 干 , 湿 季 RSBI 与 土壤 
含水 量 的 相关 性 均 高 于 NDVI, 但 在 重要 性 结果 中 ， 
NDVI(13.2%) 重 要 性 高 于 RSBI(10.0%)。 表 明 相 关 
性 较 低 的 NDVI 也 能 为 干旱 区 土壤 水 分 降 尺 度 模型 
的 构建 提供 一 些 信息 ,但 改进 后 得 到 的 EMSAVI 相 
较 于 NDVI 更 适合 作为 干旱 区 土壤 水 分 降 尺度 因子 
构建 模型 。 并 且 ,重要 性 结果 对 环境 因子 与 土壤 水 
分 关系 的 逆向 推理 具有 重要 作用 。 

(3) LST NDVI FI DEM 作为 传统 土壤 水 分 
降 尺 度 因子 ,同时 在 此 基础 上 增加 RSBI 与 引入 
SWIR 波段 改进 得 到 的 EMSAVI 构 建 干旱 区 土壤 水 
分 降 尺度 模型 。 新 旧 模 型 总 体 精度 验证 结果 表明 ， 
新 构建 的 模型 精度 尼 较 传统 模型 提升 了 0.031 ,达到 
0.913; RMSE 下 降 了 0.0013 cm .cm ,达到 0.0069 
cm tcms 

(4) 空间 一 致 性 评价 结果 表明 ,干旱 区 土壤 水 
分 降 尺 度 模型 (LBP.,=0.0585 ) 性 能 优 于 传统 降 尺 度 
模型 (LBP。 ,=0.0645 ) ,能 够 得 到 与 SMAP 产 品 之 间 
的 空间 分 布 特征 更 为 一 致 的 降 尺 度数 据 。 

本 研究 在 传统 降 尺 度 模 型 基础 上 , D n 
植被 和 裸 沙 的 环境 特征 具有 较 好 关联 性 的 光谱 指 
数 ,得 到 精度 及 空间 一 致 性 较 高 的 干旱 区 高 时 空 分 
辨 率 土壤 水 分 数据 ,能 够 对 干旱 区 土壤 水 分 降 尺 度 
研究 提供 思路 ,为 荒漠 生态 系统 碳 汇 研究 提供 重要 
的 数据 支撑 。 
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Downscaling analysis of SMAP soil moisture products in Gurbantunggut Desert 
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(1. College of Resources and Environment, Xinjiang Agricultural University, Urumqi 830052, Xinjiang, China; 
2. Xinjiang Soil and Plant Ecological Process Laboratory, Urumqi 830052, Xinjiang, China) 


Abstract: The low spatial resolution of SMAP products limits its applicability to sparsely vegetated arid regions 
and deserts with high surface heterogeneity. Considering the special environmental characteristics of sparsely 
vegetated desert areas in arid regions, traditional downscaling methods such as land surface temperature (LST), 
normalized difference vegetation index (NDVI), and digital elevation model (DEM) have been used, among 
others. Based on the scale factor, the enhanced modified soil-adjusted vegetation index (Enhanced Modified Soil- 
Adjusted Vegetation Index, EMSAVI) and the ratio sand brightness index (RSBI), which are more correlated with 
the desert surface soil moisture, were added to reflect the study area. For the downscaling factors of vegetation 
coverage and bare sand distribution, the random forest (RF) algorithm was used to build a soil moisture 
downscaling model in arid areas. The results showed the following: (1) Correlation analysis showed that 
EMSAVI (ra, ^ —0.37, Fwi = —0.34) and RSBI (ra, ^ —0.42, rw — —0.25) were good indicators of desert soil 
moisture, being superior to NDVI (ra, = -0.21, rw = 0.08). (2) The importance of EMSAVI and NDVI was 18.7% 
and 13.2%, respectively, and EMSAVI contributed more to the construction of the downscaling model. (3) The 
results obtained from the soil moisture downscaling model in dry and wet season arid regions and R of the SMAP 
product reached 0.916 and 0.910, and the RMSE reached 0.0075 cm - cm ? and 0.0063 cm + cm”, respectively, which 
are lower than the RMSE of the traditional model of 0.0013 cm°- cm”. (4) By calculating the difference (LBPc) of 
LBP (local binary patterns) to evaluate the spatial consistency, the result of the newly constructed downscaling 
model (0.0585) was better than that of traditional downscaling (0.0645). This research shows that introduction of 
the short- wave infrared band into the EMSAVI established by the vegetation index enables its better application 
to the study of soil moisture downscaling in sparsely vegetated desert areas in arid regions. 
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